Robustness of Fractional Quantum Hall States with Dipolar Atoms in Artificial Gauge 
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The robustness of fractional quantum Hall states is measured as the energy gap separating the 
Laughlin ground-state from excitations. Using thermodynamic approximations for the correlation 
functions of the Laughlin state and the quasihole state, we evaluate the gap in a two-dimensional 
system of dipolar atoms exposed to an artificial gauge field. For Abelian fields, our results agree 
well with the results of exact diagonalization for small systems, but indicate that the large value 
of the gap predicted in [Phys. Rev. Lett. 94, 070404 (2005)] was overestimated. However, we are 
able to show that the small gap found in the Abelian scenario is dramatically increased if we turn 
to non-Abelian fields squeezing the Landau levels. 



I. INTRODUCTION 

Thirty years after its discovery in GaAs-AlGaAs het- 
erojunctions [1], the fractional quantum Hall effect 
(FQHE) still remains to be the topic of current research, 
as nowadays this former sohd-state phenomenon is re- 
discovered within the scope of quantum gases in two di- 
mensions. Proposals to achieve the crucial ingredient, 
a perpendicular magnetic field, within such systems of 
neutral particles are at hand: Conceptually the sim- 
plest of them is a rotation of the system 0, |^, where 
the strength of the artificial magnetic field can be tuned 
by the frequency of rotation. However, addressing the 
regime where only the lowest Landau level (LLL) is oc- 
cupied, and at the same time guaranteeing stability of the 
atomic cloud requires a delicate balance between rotation 
and trap frequency, making an experimental realization 
of this proposal extremely hard. 

An alternative way which circumvents this problem 
is to generate artificial fields by implementing a laser- 
assisted tunneling scheme within an optical lattice '3]. 
The idea is that during a hopping process stimulated by a 
laser field, the atom picks up the phase of the laser, which 
effectively simulates the action of a magnetic vector po- 
tential. In a very similar way, a J7(l) Berry phase mim- 
icking a gauge field can be inscribed into continuous sys- 
tems via atom-laser coupling Furthermore, if the cou- 
pling involves more than two atomic states, it is possible 
to create a space-dependent degenerate subspace, which 
can be understood in terms of a non-Abelian Berry phase 
communicating between the degenerate atomic states. In 
this way, the atom-laser coupling scheme is generalized 
for synthesizing also non-Abelian gauge fields [g]. The 
same can be achieved in optical lattices, if the laser- 
assisted tunneling is made sensitive on an additional, in- 
ternal degree of freedom of the atoms 7]. These new 
possibilities have directed the attention to FQHE in such 
non-Abelian scenarios [1] . The practical feasibility of us- 



ing lasers to implement artificial magnetic fields has been 
shown in pioneering experiments [9|. Furthermore, very 
recently the implementation of spin-orbit coupling within 
a Bose-Einstein condensate can be considered an SU{2) 
gauge field realized in cold atoms (lOj . 

The main motivation for seeking new realizations of 
the old fractional quantum Hall effect is given by the in- 
triguing quasiparticles which occur as the excitations of 
fractional quantum Hall systems: Being neither bosons 
nor fermions, these so-called anyons behave exotically un- 
der interchange as instead of gaining simply a sign 
factor, a phase is obtained, which in degenerate anyonic 
systems may even be an element of some non-Abelian 
group. Due to this property, together with their topo- 
logical and thus robust nature, quantum Hall states are 
especially interesting for quantum computation [T2| . 

It is therefore most relevant to find fractional quantum 
Hall systems, where a large energy gap separating the 
ground state from the excited states guarantees a high 
degree of robustness. 

Considering dipolar atoms in a rotating trap, Ref. Q 
claims to have achieved this. In this Paper, we consider 
once again this scenario and find a much smaller gap. 
This finding directs our attention to the question if a ro- 
bust FQHE might instead be realizable in a non-Abelian 
environment. From Refs. [isj it is known that in non- 
interacting systems the non-Abelian gauge field strongly 
influences the nature of the integer quantum Hall effect. 
By choosing an anisotropic field configuration it becomes 
possible to squeeze the Landau levels. Applying to such 
a squeezed scenario the thermodynamic approximation 
14 j used in the Abelian case in Ref. Q , we show in this 



Paper, that a dramatic increase of the energy gap can be 
achieved. 

As in dipolar systems it turns out that the gap cru- 
cially depend on how it is defined, we first discuss this 
matter in Section [Hi which we conclude with a precise 
definition of the gap. Subsequently, we describe in Sec- 
tion mil an analytic approximation allowing to evaluate 
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this gap in the thermodynamic hmit. Section [IVI is ded- 
icated to specifying the gauge potential. Here we show 
that Abehan and non-AbeUan fields can be treated on 
the same footing. Finally, the concrete calculation and 
its results are described in Section IVl 



II. FRACTIONAL QUANTUM HALL EFFECT 
AND DEFINITION OF THE GAP 

The physics of non-interacting particles confined to a 
plane under the influence of a perpendicular magnetic 
field is understood in terms of the quantized Landau lev- 
els. While this scenario may give rise to an integer quan- 
tum Hall effect, the FQHE requires repulsive interactions 
between the particles as its crucial ingredient. In the 
context of FQHE as a solid-state phenomenon, almost 
exclusively Coulomb interactions have been considered. 
For neutral fermions, however, it is most likely to have 
dipolar interactions. We therefore study the following 
Hamiltonian (h ~ 1): 

N ^ N 

H = E 2^;^ [P^- - ^(^^O]' + E ^dd(r-„ r,). (1) 

j = l j<k 

Here m is the mass of the particles and Vdd{rj,rk) — 
cP/\rj — Tfcl^ the dipolar interaction between the parti- 
cles with dipolar moment d. The gauge potential A is 
supposed to describe a constant gauge field perpendicu- 
lar to the two-dimensional system, and will be specified 
in Section HVl 

It is known that for practically any repulsive inter- 
action a very good trial wave function for the ground- 
state is the so-called Laughlin function J^], which can 
be found by filling the lowest Landau level up to a filling 
factor v 

N / ^ \ 

*l({z,}) = Mllizk - z,fl^ exp - ^ I Ml ) , 

(2) 

where zi = Xi + iyi are the positions of the particles, Af is 
a normalization factor, and Iq = y^TjBo is the magnetic 
length corresponding to a magnetic field of strength Bq. 
If 1 /u is odd, the function is fully antisymmetric and thus 
describes fermionic systems, while the opposite is true for 
even l/v. 

The Laughlin ground-state can be considered a homo- 
geneous liquid state. It is excited by piercing a hole into 
it. Choosing for simplicity the quasihole to be at the 



origin, the excited state is described by 

N 

*qh({^,})=AA'n^J*L, (3) 

with A/"' a normalization factor. 

A gap, vaguely defined as the energy difference sepa- 
rating the Laughlin state from the quasihole state, guar- 
antees the robustness of the FQHE. But as one may think 
of different ways of creating a quasihole, it is important 
to be precise in the definition. In a finite system, a quasi- 
hole can be created in three different ways [l3l • 

(i) by reducing the particle number, 

(ii) by increasing the area of the system at constant 
particle number, 

(iii) by changing the magnetic field at constant particle 
number and constant area. 

In the electronic system, in addition to the Coulomb re- 
pulsion between the electrons, a positively charged back- 
ground (originated from ions) is present as a consequence 
of the electroneutrality of the system. The corresponding 
background potential stabilizes the Coulomb system, but 
it also adds to the energy of the quasihole, compensat- 
ing the losses in the energy of the direct electron-electron 
Coulomb interaction due to lowering the number of elec- 
trons or their density. As a result, for any choice of the 
definition a positive energy gap is obtained [13, ■ 

Contrary to the electron case, dipolar systems have 
no such background potential, and the definition of the 
gap is crucial. As we will explicitly see in Section [Vj 
energy "gaps" defined according to (i) or (ii) have a neg- 
ative sign. This would mean that the system is unstable 
against the creation of quasiholes. In fact, this, at first 
view, counter-intuitive result is not astonishing: Since (i) 
and (ii) results in diluting the system and, therefore, in 
lowering the potential energy, it is obvious that in these 
cases we should find an energy gain. Therefore, the only 
meaningful definition for the energy gap in a dipolar sys- 
tem is according to (iii) , where we compare the energy of 
a quasihole with the ground state energy of the system 
with the same number of particles occupying the same 
area - similar to the definition of the quasihole excitation 
energy in a normal Fermi system. 

III. THERMODYNAMIC APPROXIMATION 

In order to evaluate the gap, we follow an approach 
developed in Refs. [ij. Based on the plasma analogy 
, relating the physics of the Laughlin state to the one 
of a classical one-component plasma, analytic expressions 
for the correlation functions gQ{zi,Z2) and gqh(zi,22) of 
the Laughlin state and the quasihole state can be derived: 
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Note that in these expressions, the quasihole state is re- 
lated to the Laughhn state by reducing the particle den- 
sity in the center of a thermodynamically large system, 
which according to the classification of the previous sec- 
tion is a quasihole of type (i). 

Given these expressions, we may evaluate the energy 
difference A between these two states by solving the fol- 
lowing integral 



[9qh{zi,Z2) - go{zi,Z2)] ■ 



(7) 



Fortunately, we are able to derive a relation between 
the quantity from Eq. ^ and the gap related to a 
quasihole state according to (iii). First, we notice that 
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Eq^\ where these are the energies of 



N — V particles in a quasihole state, and N particles in 
the Laughlin state. Now we may write Eq — iVeo, where 
eo is the energy of one particle in the Laughlin state. 
This quantity can readily be evaluated by substituting 
z- = zi — Z2 in Vdd(zi, Z2) = V{z-) and z_ = zi — Z2 in 
go{zi,Z2) = go{z-), and integrating: 



2v 



dz_ gQ{z_)Vdd{z-). 



(8) 



Demanding a constant particle number, we may re- 
define the gap as 
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qh 



-E, 



{N-u) _ 
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(9) 



This definition describes a quasihole created according 
to (ii). However, as argued in the previous section, apart 
from a fixed particle number, we should also demand 
a fixed volume. Therefore we notice that each Landau 
state occupies an area a — ttIq. For the Laughlin state 
with N particles to occupy the same area as the quasihole 
state with N particles and one quasihole, we thus have 
to modify the magnetic length I'q of the excited state 
according to: 
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(10) 



Now we have to note that the energies in the dipolar 
system scale with 1^^ . Since we wish to compare states 
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at different magnetic fields, we define the gap at constant 
particle number and constant volume as: 
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By noticing that 



TV 



(11) 



(12) 



and approximating N/ {N — zy) w {N + v) /N for large N , 
we find with Eq. pHI) : 



Ay = A + -j/eo- 



(13) 



Before we calculate this quantity, we discuss different 
gauge potentials that might be realized. 



IV. 



DIFFERENT GAUGE POTENTIALS 



For the discussion of FQHE, it is feasible to choose a 
symmetric gauge for the gauge field. A perpendicular, 
magnetic field of strength Bq is then described by the 
magnetic vector potential ^mag('''j) = BQ{—y,x). In the 
context of artificial fields, it is possible to generalize this 
potential to be an element of SU{2). We choose it to 
have the form: 



A{rj) = So(-y,x) + {aay,(iax), 



(14) 



where a, and /? are additional, controllable parameters, 
and Ox^y are Pauli matrices. Proposals to realize such 
auge potentials have been made both for lattice systems 
and for trapped gases coupled to laser fields Q. Eq. 
contains the limit of a magnetic field, as can be seen 
by choosing a — (3 ^i). For finite a and /3, the potential 
([T^ yields a constant non-Abelian gauge field perpendic- 
ular to the system. For the isotropic configuration a — (3 
it has been shown in Ref. Q that this potential splits 
the Landau levels, which, by assuming a short-range re- 
pulsion between the particles, may give rise to a FQHE 
with non-Abelian anyons. The most general configura- 
tion with a ^ P has been considered for non-interacting 
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lattice systems in Ref. [Tj]- Then the SU{2) gauge po- 
tential produces an anisotropic space-time. Giving rise 
to Dirac points in the bandstructure of the free system, 
we may measure this anisotropy as the ratio of the sound 
velocities in x- and y-direction: Cx/cy = \a\/\/3\. In Ref. 
[l3l | it has been shown that this situation is best described 
by introducing a squeezing parameter 



-tanh 



— If y ^ 



(15) 



and replacing the original variable z by squeezed one z: 



iy] [z{£,) = cosh^ z — sinhf z] . 



(16) 



We are then able to treat both the Abelian scenario and 
the non- Abelian scenario on the same footing: The gener- 
alization of the Laughlin wavefunction and the quasihole 
wavefunction, Eqs. © and ([3]), to systems with a non- 
Abelian field are straightforwardly given by making the 
replacement Eq. (|16p . Accordingly, the correlation func- 
tions derived for the states Eqs. © and ([3]), also hold 
for the corresponding squeezed states, if we again make 
the substitution Eq. ([TC]). The gap, as defined in Eq. ([7]) 
can then be evaluated by the integral 

A(e) =i J d^Zi I d^Z2 Vm{zi,Z2) (17) 

[5qh(^i(o,2~2(e))-5o(5i (0,^2(0)]- 

In the same way, we generalize the ground-state energy 
defined in Eq. dH) to be a function eo(^) of the squeezing. 
Following the derivation as in Section lTTll we finally arrive 
at the equation Ay(^) = A(^) + |i^eo(0- 

V. RESULTS 

While in the Abelian scenario with ^ = parts of the 
calculation can be done analytically, the squeezed sce- 
nario is more complicated as it demands a fully numerical 
treatment. Since in both cases the steps of the calculation 
are the same, we describe in details only the procedure 
for the Abelian case. 

The main difficulty consists in evaluating the integral 
Eq. [71 First we have to specify the coefficients Cj in Eqs. 
dH) and It is shown in Ref. [ll| that by setting all 
Cj = 0, a system with a completely filled Landau level is 
described, v = 1. For this choice of Cj, the resulting cor- 
relation functions are denoted by g^^^ and g^^ . In order 
to have a FQHE, we need a fractional filling, v — 1/q, 
which requires the coefficients Cj with j < q to be non- 
zero. For fermions, the most robust effect is expected for 
1/ = 1/3, where the choice Ci = 1 and C3 = —1/2 is 
best suited. We call the corresponding correlation func- 
tions gl^^ and g^^ , and also define for convenience the 

differences Sq = 5o^^ - ffo^^ and Sqh = ffqh - .9qh • 

Turning now to the integral Eq. ([7]) with i> = 1/3, 
we note that in the Abelian limit it reduces to the one 



considered in Ref. ^. As our numerical result, however, 
drastically differs from Ref. [31, a careful analysis is of 
order. Therefore we split the integral Eq. ([7]) into two 

parts. Pi ^ J dzi J Azi V;3d(5qh ~ 5o^'')' which is ana- 
lytically solvable, and P2 = J dzi J dz2 Vdd(Sqh — Sq^-*), 
which we treat numerically. 

For the analytic part we find Pi = —^f^T^jv^y^. Note 

that for V — \, this negative number would be the full, 
completely analytic result for the energy difference A de- 
fined in Eq. ([7]) . This clearly shows what we have antici- 
pated in Section ini namely that this definition is not the 
appropriate one for the energy gap in a dipolar system. 

Before we evaluate P2 numerically, we examine the 
asymptotic behavior of the integrand. As the divergence 
in the interaction term for zi — >■ Z2 is compensated by 
the vanishing of the correlations, this limit can easily be 
handled by a regularization of the integral. The limit 
of = 2:1 -|- Z2 — >■ 00, however, turns out to be prob- 
lematic: For finite particle distance, \z\ — Z2\ < 00, this 
contribution is not suppressed by the interaction, and 
the convergence of the integral Eq. ([7]) requires that Eq 
and Sqh have the same asymptotic behavior. However, 
the completely different structure of both functions ob- 
scure the latter. Contrariwise, we should note that if we 
truncate the infinite sum in Eqh, this expression gets ex- 
ponentially damped for large center-of-mass coordinates, 
while Eq depends only on the relative coordinates, yield- 
ing So - Eqh ^ for |z+| -J> 00. 

To circumvent this problem, we bring Sq to a form sim- 
ilar as Eqh, which is possible by factoring out a damping 
exp[— (|zip -I- |z2p)/2] and Taylor expanding the remain- 
ing exponential exp[|z_|_p]. We are then able to write 

' ^ll^ + U-.l^ ^ -2Q ^ |F(°'(zi,Z2)|2 



3 ■' k=G 



(18) 



(0) 



k k 



3,k 



(19) 



r=0 s=0 



As now each term in both Eq and Eqh is damped by a 
factor exp[— (|zip -I- |z2p)/2], they all vanish in the limit 
|z+| — 00, and we may truncate the infinite sums at a 
sufficiently large value of k. Note that due to the different 
orders in zi and Z2 of Fj^k in Eq. ^ and in Eq. p^ . 
the sum in Eq should contain two more terms than the 
sum in Eqh for a quick convergence. 

Now we are able to perform the numerical integration. 
The error due to the truncation still is 5% for 10 terms, 
but can be minimized by a finite-size analysis of our re- 
sults. We then find P2 = (0.1875 ± 0.0010)^, where 
the numerical error has been approximated by the devi- 
ation from the smooth fit in Fig. [TJ With this, we find 
A = 0.5(Pi + P2) = -(0.0455 ± 0.0010)#. As we have 
argued in Section [III this negative value is due to the 
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FIG. 1. The gap Av at constant volume and constant particle 
number as a function of the squeezing parameter ^ : The dots 
are obtained by a numerical evaluation of Eq. ([7]) for different 
^, a fit of this data yields the solid line. 

reduced density of the system. 

We continue with calculating the gap as defined in 
Eq. (IH). Therefore we have to evaluate the integral 
Eq. which in the isotropic case, ^ = 0, reads 

eo = ~ -gf^ jr, from which we find that also 

Ajv < 0. The negative values for A^v can be under- 
stood by noticing that as the particle taken away at the 
origin now has been added at the edge of the system in- 
creasing its volume, so A^r corresponds to the energy of 
a quasihole according to (ii). As long as such a process 
is possible, the system is unstable as it tries to reduce its 
density by diluting. 

Finally, we turn to definition Eq. (fT3|) . Only in this 

case, we obtain a positive gap. Ay = (0.0132±0. 0020)4-, 

9 

which however is much smaller than the number found in 
Ref. [1, (0.9271 ± 0.019)|^, but compares weU with the 

'o 

gap obtained via exact diagonalization of a small dipolar 
system in Ref. [ll], where the discrepancy to Ref. [1] 
has been attributed to the different system size. 

Repeating this procedure for finite squeezing ^, we have 
to evaluate A(^) as defined in Eq. ([T7]), but now the 
whole integral has to be solved numerically. Again we 



find negative values for A(^) and Ajv(^), which even de- 
crease with larger ^. However, as also the ground-state 
energy eo(^) increases with ^, the gap Av(^) at constant 
particle number and constant volume finally has a pos- 
itive balance for all ^. As shown in Fig. [l] it increases 
with ^, and a convenient fit to the numerical data is found 
to be: 

Av(0 = Ay(0)exp(aC'). (20) 

We obtain a = 0.529 and Ay(0) = 0.0152 d^/ll 

To understand this behavior, we note that the squeez- 
ing allows the particles to get closer in one direction, 
while the particle distance is increased in the other di- 
rection. Due to 1/r^ behavior of the dipole-dipole in- 
teraction, the interaction energy is much more sensitive 
to changes of the density distribution at short distances 
rather than at large ones. Thus, compressing in one di- 
rection and stretching in another one increases the inter- 
action energy. As a consequence of Eq. p^ . this gives 
rise to a bigger energy gap. 



VI. CONCLUSION 

Summarizing this work, we have shown that dramat- 
ically differing from the predictions in Ref. 3, only 
a small energy gap separates the Laughlin state from 
quasihole excitations in systems of dipolar quantum gases 
with artificial magnetic fields. However, by considering 
scenarios where a non-Abelian gauge field introduces an 
anisotropy into the system, an exponential increase of the 
gap can be achieved, and may allow for robust fractional 
quantum Hall states. 
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